A MEAN VALUE THEOREM FOR THE 
DIOPHANTINE EQUATION axy - x - y = n 
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Abstract. In this paper, we prove an asymptotic formula for the 
average number of solutions to the Diophantine equation axy — x — 
y = n in which a is fixed and n varies. 



1. Introduction 

People has been considering Diophantine equations involving products 
and sums of some variables for a long time. The Diophantine equation 

k k 

\\ Xj - ^ Xj = n (1) 

8=1 i=l 

was studied by various people during the past a few decades. It is easy 
to see that there always exists a few trivial solutions with most of Xj's 
equal to 1. So people are asking about the number of solutions of this 
equation with all x^ > 1. 

The case when n = is very special, since it concerns the number of 
/c-tuples with equal sum and product. In this case, it is conjectured by 
Misiurewicz [2] that k — 2,3, 4, 6, 24, 114, 174 and 444 are the only val- 
ues of k for which there are only trivial solutions. For general n, very 
little is known except that in 1970s Viola [6] proved that if Ek(N) de- 
notes the number of positive integers n < N for which ([Q) is not soluble 
in integers Xi, x 2 , ■ ■ ■ , > 1 then E k (N) = iVexp(— Cfc(log jV) 1-1 /( fc+1 )) 
for some positive constant Cfc. It is believed that for large n equation (JT]) 
always has a nontrivial solution, which nevertheless is an open question 
in this area. 

On the other hand, the case that k = 3 has received extensive at- 
tention, and several variations of this problem were studied. Brian 
Conrey asked whether the number of solutions in positive integers to 
the equation 

xyz + x + y = n (2) 
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can be bounded by £ (n e ) for any e > 0. Kevin Ford posed a general- 
isation of this problem, in which one would like to show that there are 
E (\AB\ £ ) nontrivial positive integer solutions to the equation 

xyz = A(x +y) + B (3) 

for given nonzero A, B G Z. 

In this paper, we consider another variation of the case that k = 2, 
namely the following equation 

axy — x — y = n (4) 

where a is a positive integer and n is any nonnegative integer. This can 
be viewed as equation ([3]) in which z is fixed and A = 1. Hence if the 
number of solutions of equation (j3]) is well understood, then one can 
probably understand the number of solutions of equation ([3]) simply by 
averaging over a. 
Let 

R a {n) = Card {(x, y) 6 N 2 : axy — x — y = n} . 

Here we are considering the number of positive integer solutions of 
equation (j2J) when a is fixed and n varies. A sharp asymptotic formula 
is established in this paper on the average of R a (n) over n. Notice the 
case that a = 1 is trivial, since then Ri(n) = d(n + l) is just the divisor 
function of n + 1, the average of which is relatively well understood. 

Theorem 1. For positive integers a > 1 and N > 1, we have 
Ra(n) = ^(N log N - C(a)N^) + A a (N) 

0<n<N 

where 

C(a) = 2^M + 2^^ + loga + 2 7 + l (5) 

^ a ' p\a ^ 

and 

A a (iV) « 0(a)^(log(aiV)) 2 . (6) 

Here = / + °° e~ t t s ~ 1 dt is the standard T function, and 7 is the 
Euler constant. 

In fact, since the error term above is roughly of size V aN ^ log(aiV) j , 

it is conceivable that the main term will be inferior to the error term 
when a ^> N3, So in order for the above asymptotic formula to really 
make sense, one would impose a condition on a, such as a iV3 / log N. 

Moreover, one can argue what is the right order of magnitude of the 
error A a (N). In view of R\{n) = d{n + 1), one can think R a {n) as a 
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"generalized" divisor function. Hence Theorem [T] just proves a mean 
value theorem for such a "generalized" divisor function. Since for the 
classical divisor function, the error is believed to be 0{N l / A+£ ). It is 
very natural to pose such a conjecture for our error A a (N). The author 
suspects that following the van der Corput method on exponential sums 
as in the classical case, one can show A a (iV) = O a (A rl//3 ~' 5 ) for some 
5 > 0. 

Remark. It's not hard to adapt the method in this paper in order to 
deal with equations like 

axy — bx — cy = n 
and prove similar asymptotic formulas. 

2. Preliminary Lemmas 

We state several lemmas before embarking on the proof of Theorem 
1. The content of Lemma 2 can be found, for example, in Corollary 1.17 
and Theorem 6.7 of Montgomery & Vaughan j3], and Lemma 3 can be 
deduced from Theorem 4.15 of Titchmarsh [5] with x = y = (\t\/2n) 1 ^ 2 . 

Lemma 1. When a > 1 and \t\ > 2, we have 

—^—7 < C{a + it) < log 
log |i| 

Lemma 2. When < a < 1 and \t\ > 2, we have 

((<t + it) < |t|^log(|t|). 

Lemma 3. Let \ be a non-principle character modulo a and s = cr + it 
and assume that tel. Then 

L(s, x) < log(a(2 + when a > 1 

and 

L(s, x) < (altl) 1 ^^ , when - < a < 1. 

Proof. The first part follows from Lemma 10.15 of MV [4] . Now suppose 
that x is primitive. Then by Corollary 10.10 of MV jl], 

L( S , X )«(a|t|)^ <T log(a(2+|t|)) 

when o~ < 0. Then by the convexity principle for Dirichlet series, for 
example as described in Titchmarsh [5] (cf. exercise 10.1.19 of MV [1]), 

L( S , X )«(a|t|) i ^ +£ 
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when < o < 1. The proof is completed by observing that if \ < a < 1 
and x modulo a is induced by the primitive character x* with conductor 
q, then 



L(s, x) = L(s, x*) U(l ~ X*(p)p- S ) « \L(s, X *)|2 



uj(a) 



p\a 



Lemma 4. Let T > 2, then we have 



□ 



£ / \L(Uit, X )\*dt«^-TlogT. 

mod a 

A proof of this lemma can be found for example in Montgomery [3] . 

Lemma 5. Let a be a positive integer greater than 1 and w > 0, we 
have 



- = - logw - - logo ) + 0(l/w) 

n<w v v a > 

n=—l mod a 



E 



Proof. By Abel summation, the left hand side above is 

E 1 = 



n<w 
n=— 1 mod a 



10 + 1 


1 

- 4 




t+1 


a 


w 


Jl 


a 



1 



t + 1 



a 7 X at 2 
' log u> 4 2 



£ + 1 1 tft 



a J t 2 
t + 1 1 dt 



a t 2 



+ o(ih 

+ 0(l/tw) 



Recall that the digamma function is defined ), and ip'(z) 



^00 1 



has a series expansion E^ fc=0 ( z +k) 2 



. So 



HI 



w-i)f = E 

h=0 
1 



r 4 1 
a 

r 4 1 



- W - 1 



(ah 4 ?") 2 



dr 
(7) 
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Notice that 



r + 1 



( 0, if 
1, if 



< r < 1 

1 < r < 2 



{r} - 1 = < 



a — 2, if a — 2 < r < a — 1 
—1, if a — 1 < r < a 



Hence ([7j) is equal to 

'o-2 „i 



( a ) a^^(a) 



, /a — 1 + r 



V> - 
a 



^(1)-V 



dr 
a — 1 



a — 1 



+ log a + 7 



The last equality follows from a well known property of the digamma 
function ip. Now the lemma is established after the observation 7 = 

2 _ iii+irft. □ 

Lemma 6. Lei a be a positive integer greater than 1, then we have 
1 



6(a) ' vv / v "w a \ T 

Xt^XO 
mod a 



1 fH*?), ~bgp, . , 



Proof. Let u> be large compared to a (eventually we will let w goes to 
00). Then for non-principal characters x modulo a, by Abel summation 



£(i,x) = E 



X[n) 



Hence 



n<w 



W) £x(-i«i,x) 

vv ; X^Xo 

mod a 



0(a/w). 



0(a) ^— ' n 

Xt^XO ™<™ 
mod a 
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The main term on the right is 

— y x(-i)y^--— y -■ 



6(a) ^— ' ^— ' n 0(a) ^— ' n 

mod a (n,a)=l 



We have 



— — - ( log( w/m) + 7 + Oirn/w) 
( log W + ^^- + 7 ) +0(d(a)/ W ). 

V p\a P J 7 



n ^ — ' m z — ' n 

n<w m\a n<w/m 

(n,a)=l 



m 

m\a 



6(a) 



p\c 

Here we are using the fact that — Y] m \„ ^£r- logm = Y]„,„ this 

o z-^m\a m ° a t-^p\a p—1 ' 

is because 

m|a p|a fe|a/p 

(p,fc)=l 



p\a p \a 

v'+v 

s ( i - i J n ( 1 ~ ^ 

p\a y \ Vj p'\a V y 

0(a) ^ logp 



a ^ p - 1 

p|a 



On the other hand, we have 

jzs y xi-^y = y - 



6(a) ^— ' ^— ' n ^— ' n 

mod a iiE-lmoda 



And by lemma [5], this is 

1 (, Y ,(a - X 
a 



~ [ logw- - log a ) t()(1/ic). 
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Thus we have shown that 

^ X^O 
mod a 



+ 7 +0(a/w) 



-(log«, - ^^ry - loga) - - I logw + E^TT 

V a ; \ p|a y 

^ a ' p\a J 

Now the lemma is established when we let w — > oo in the above. □ 

3. Proof of Theorem Q] 

The starting point of the proof is the following observation. One can 
rewrite equation (j3J) in the following form 

(ax — l)(ay — 1) = an + 1. (8) 

Namely we are going to count the following quantities, 

R a (n) = Card {(x,y) G N 2 : (ax - l)(ay - 1) = an + l} 

and 

S (iV) = £ R a (n). 

0<n<N 

After the change of variables u = ax — 1 and i> = ay — 1, it follows 
that R a (n) is the number of ordered pairs of natural numbers u, v such 
that uv = an + 1 and w = t> = — 1 (mod a). 

Now the residue class u = — 1 (mod a) and t> = — 1 (mod a) are 
readily isolated wa the orthogonality of the Dirichlet characters x mod- 
ulo a. Thus we have 

= £ £ 1 

0<n<N uv=an+l 
us— 1 mod a 
u= — 1 mod a 

= £ £ i 

m<M ™=m 
— ms— 1 mod a 



d=— 1 mod a 



^ £ £ *<-U*(-i) £ £ 

T v ' Xi 



Xi X2 m<M uv=m 

mod a mod a 



where M = aiV + 1. 
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Let 



Then we have 



^ ^ Xi X2 m<M 
mod a mod a 

We analyze this expression through the properties of the Dirichlet series 
fxuxM = £ a? " (Xl ; X2) = L( S ,Xi)^,X 2 )- (9) 

77. 

m=l 

This affords an analytic continuation of f XuX2 to the whole complex 
plane. 

By a quantitative version of Perron's formula, as in Theorem 5.2 of 
MV [1] for example, we obtain 

1 f ao+iT M s 

a m{Xl, X2) = 7T- / /xi,X 2 ( S ) ds + R {Xl, X2), 

where do > 1 and 

R{Xl,X2) < Ef <m<2M |Om(Xl»X2)| mm (l, ^ M 

I 4fg+Mfj v^oo |q m (xi»X2)| 

Here means that when M is an integer, the term clm{Xi>X2) is 
counted with weight |. 

Let 0- o = 1 + j^g. By © we have |a m (xi,X2)| < d(n). Thus 

E^<<CW<<(logB) 2 

m=l 

and so -R(xi? X2) M 1+e T _1 , for any e > 0. Hence 

E Mxi, ^) = T~i fxuxM—ds + 0(( — + l\M s \. 

The error term here is 

< M £ 

provided that 

T > M. 
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The integrand is a meromorphic function in the complex plane and is 
analytic for all s with > | except for a possible pole of finite order 
at s = 1. Suppose that T > 4. By the residue theorem we have 

1 r +lT t < j 



2ni 



Jao-iT J\-iT J\+iT 



L(s,xi)L(s, X 2)M s ^ 



2 

+ Res s=1 ( L(s,Xi)L(s,X2)— j ■ 

Hence, by Lemmas 1, 2 and 3, the contribution from the horizontal 
paths is 

< T-^aTfM + T^^aTf'^M 1 ' 2 
and provided that T > M 5 this is 

< M-\ 

On the other hand, the contribution from the vertical path on the 
right is bounded by 

M^2 4 / \L(l + it,Xi)L(l+it, X 2)\dt. 

J 2 k 

2 k <T 

And by Lemma 4 

f u - f u 1 ( l+iT L{s,Xi)L{s,X2)M° J 
^ Xi(-l)x 2 (-l) — ds 

XUX2 Am H-iT S 

mod a 

2 k+1 / \ 2 

dt 



2 fe< T ^ 2fc V X / 
mod a 

r 2k+1 

«M^2" fc 0(a)^ / |L(| + zt, X )| 2 rft 



2 fe <T ^ 

mod a 

^3/ 



2 k <T 

6 3 (a) i . 
< I-^M2(logM) 



a 
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on taking 

T = M 5 . 

Hence we obtain 



1 / M s \ 

Sa{N) = ^{a) E E Xi(-l)X2(-l)Res s=1 ^ Xl , X2 ( S )— J+A a (iV) 



XI X2 
mod a mod a 



where 



A a (iV) < ^^(log Mf < 4>{a)J^(\og{aN)} \ (10) 

It remains to compute the residue at s = 1. 
By (JH]) there are naturally two cases, namely 

(i) xi = X2 = xo; 

(ii) only one of Xi an d Xi is equal to xo while the other one is equal 
to x Xo- 

In the latter case the integrand has a simple pole at s — 1 and the 
residue is 



J] (l - -) L(l, X )(aN + 1) = <P(a)L(l, X )N + 



x). 



By lemma El the sum over x f° r the second term above is small, hence 
can be absorbed in A a (N). While in the former case, the integrand 
has a double pole at s — 1 and the residue is 

2 



p\a 

Hence we have shown that 



log M — M 



S a (N)= ±((aN+l)\og(aN+l)-aN-l) 
+ (l^E x(-mi,x))N + A a (N). 



Xt^XO 
mod a 



Now by lemma El this is 

- (N log N - C(a)N) + A a (N) 

where C(a) and A a (N) are given by (JSJ) and dfJJ) respectively. 
This completes the proof of Theorem [TJ 

□ 
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